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ABSTRACT 

Context. The study of the Luminosity Function (LF) of Lyman Break Galaxies (LBGs) at z~7 is very important for 
ascertaining their role in the reionization of the Universe. These galaxies can be used also to investigate in detail the 
processes of formation and evolution of galactic structures in the infancy of our Universe. 

Aims. In this work we plan to perform a detailed and critical analysis of the statistical and systematic errors in the 
2 ~ 7 LF determination. 

Methods. To this aim, we have assembled a large sample of candidate LBGs at z ~ 7 from different surveys, spanning 
a large variety of areas and depths. In particular, we have combined data from the deep (J < 27.4) and ultradeep 
(J < 29.2) surveys recently acquired with the new WFC3 NIR camera on HST, over the GOODS-ERS (~40 sq. 
arcmin.) and the HUDF (~4 sq. arcmin.) fields, with ground based surveys in wide and shallow areas from Hawk- 
I@VLT and HyperSuprimecam@Subaru. We have used public ACS images in the z band to select z-dropout galaxies, 
and other public data both in the blue (BVI) and in the red bands to reject possible low-redshift interlopers. We 
have compared our results with extensive Monte Carlo simulations to quantify the observational effects of our selection 
criteria as well as the effects of photometric scatter, color selections or the morphology of the candidates. 
Results. We have found that the number density of faint LBGs at z ~ 7 is only marginally sensitive to the color selection 
adopted, but it is strongly dependent from the assumption made on the half light distributions of the simulated galaxies, 
used to correct the observed sample for incompleteness. The slope of the faint end of the LBGs LF has thus a rather 
large uncertainty, due to the unknown distribution of physical sizes of the z ~ 7 LBGs. The implications of these 
uncertainties have been neglected by previous works. 

Conclusions. We conclude that galaxies at z ~ 7 are unable to reionize the Universe unless there is a significant evolution 
in the dumpiness of the IGM or in the escape fraction of ionising photons or, alternatively, there is a large population 
of z ~ 7 LBGs with large physical dimensions but still not detected by the present observations. 

Key words. Galaxies:distances and redshift - Galaxies: evolution - Galaxies: high redshift - Galaxies: luminosity function 
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significant variations in the escape fraction of ionising 
photons. Quantifying the number density of z > 7 galaxies 
is therefore critical in order to check whether there are 
sources responsible for the re-ionization, such as Pop III 
dominated primordial galaxies, mini-black holes or others 
(jVenkatesan et al. 20031 IMadau et al. 2004j) . 

The different results obtained for the LF of LBGs 
at z=7 (ClOa, C10b,[Ouchi et al. 200j|Capak et al. 2011 
IQesch et al. 20101 IMcLure et al. 2"0T0l IBunker et al. 2010 
Wilkins et al. 2010, Hickey et al. 2010| lYan et al. 2010 



Finkelstein et al. 2010[) should be ascribed both to the 
different data-sets used and to the different techniques 
adopted to analyse the data. We want to study in detail 
the possible systematics/uncertainties in the LF estimate, 
in particular due to the choice of the color cut adopted, 
the morphology of the z=7 galaxies assumed during sim- 
ulations, and the statistical tools used in the LF deriva- 
tion. Moreover, extending the LF at z=7 down to faint 
magnitude limits is fundamental to break current degen- 
eracies between M* and $*, and to put strong constraints 
on the number density of faint LBGs at high-z. Here we 
combine a re-analysis of the extremely deep WFC3 ob- 
servations ([Oesch et al. 20151 with the selection of bona- 
fide z ~ 7 LBGs in wide areas of the sky (ClOa, ClOb, 
IQuchi et al. 20 09) in order to derive stringent constraint to 
the reionization process at z ~ 7. 

The paper is organised as follows. Section 2 and 3 de- 
scribe the imaging data used, the photometric analysis, and 
the selection of the z ~ 7 LBG candidates. The simula- 
tions carried out to derive incompleteness and contamina- 
tion from lower-z galaxies are discussed in Section 4. In 
Section 5 we discuss the systematics affecting the LF (color 
criteria and morphology of the galaxies adopted during the 
simulations) , with our results on the LF derivation. Section 
6 and 7 are devoted to our discussions and conclusions. 

Throughout the whole paper, observed and rest-frame 
magnitudes are in the AB system, and we adopt the A- 
CDM concordance model (Hq = 70km/s/Mpc, VLm — 0.3 
and n A = 0.7). 



2. Data and Photometry 

We have analysed two different data-sets observed with 
the new NIR camera of HST, the Wide Field Camera 3 
(WFC33): the Early Release Science on the GOODS-S field 
(GOODS-ERS) and the Hubble Ultra Deep Field (HUDF) 
programs. 

2.1. ERS 

The GOODS-ERS WFC3/IR dataset is a total of 60 HST 
orbits consisting of 10 contiguous pointings in the GOODS- 
South field (HST Program ID 11359), using 3 filters per 
visit (1o98j J125, #160)7 and 2 orbits per filter (for a to- 
tal of 4800-5400s per pointing and filter). The total area 
covered by the GOODS-ERS is ~ AOsq.arcmin. down to 
Y=27.3, J=27.4, and H=27.4 magnitudes at 5er in an area 
of ~ 0.11 arcsec 2 (corresponding to 2 times the FWHM 
of the images). The ERS data were acquired and used 
to demonstrate the WFC3 capabilities on the study of 
low-z and high-z galaxies (Windhor st et al 201 ip . and the 
raw images were publicly released by the instrument team. 



http:/ /www. stsci.edu/hst/wfc3 



Thanks to its fine pixel sampling and high sensitivity 
in the NIR, WFC3 is an ideal instrument to study the 
z > 7 galaxy populations. The raw images of the ERS 
dataset have been reduced and used to search for z drop 
candidates by IWilkins et al. 20101 IWilkins et al. 20111 and 
IBouwens et al. 20111 

The WFC3 ERS images were reduced using the 
Multidrizzlc software ([Koekemoer et al. 2002|) . adopting 
on-orbit SPARS100 darks and also correcting for the gain 
differences between different quadrants. Large-scale resid- 
uals on each exposure which are apparently due to scat- 
tered light have been removed, and the satellite trails 
have been masked on affected exposures. The astrome- 
try of the images is bound to the GOODS-S v2.0 z-band 
([Giav alisco et al. 2004). The matching is robust even for a 
single exposure, with 400-600 sources matching, so that the 
overall alignment seems to be robust to a level of ~ lOmas. 
The pixel size of the reduced images is 0.06 arcsec (a factor 
of 2 rebinning was applied with respect to the original pixel 
size of 0.135x0.121 arcsec 2 ), with a PSF of 0.18 arcsec. 

2.2. HUDF 

The HUDF WFC3/IR dataset (HST Program ID 11563) 
is a total of 60 HST orbits in a single pointing 
(jOesch et al. 20101 IBouwens et al. 2010a|) in three broad- 
band filters (16 orbits in Y105, 16 in J125, and 28 in Hieo)- 
It is the deepest NIR image ever taken, reaching Y=29.3, 
J=29.2, and H=29.2 total magnitudes for point like sources 
at 5 sigma (this S/N ratio is computed in an aperture of 
~ 0.11 arcsec 2 , corresponding to one FWHM). The area 
covered by the WFC3-HUDF imaging is 4.7 sq.arcmin., 
and the IR data have been drizzled to the ACS-HUDF data 
([Beckwith et al. 2006[) . with a resulting pixel scale of 0.03 
arcsec, with a FWHM of 0.18 arcsec. 

The HUDF images were reduced in the same way as 
the ERS dataset, i.e. background effects are properly cor- 
rected (on-orbit darks, quadrant gain corrections, multi- 
plicative sky-flats and additive low-level residuals on large 
scales all removed), and the images have been astrometri- 
cally matched using the ACS HUDF sources. For the Y105 
filter, we have selected a subset of 7 visits out of 9 that 
are not badly affected by persistence from bright targets 
observed before these observations; the images in the other 
two filters (J and H) are not affected by this persistence 
problem. The magnitude limit of the cleaned Y-band image 
is 29.0, at 5 sigma in - 0.11 arcsec 2 . The GOODS-ERS and 
HUDF datasets will be described in detail in Koekemoer et 
al. (2011, in prep.). 

2.3. Photometry 

The photometric catalog of the ERS and HUDF fields 
have been derived in a consistent way. Galaxies have been 
detected in the J 12 5 band, and their total magnitudes 
have been computed using the MAG_BEST of SExtractor 
(jBertin fc Arnouts 1996|) . for galaxies more extended than 
0.11 arcsec 2 , while circular aperture photometry (diameter 
equal to 2 times the WFC3 FWHM) is used for smaller 
sources. The magnitudes which are computed in a circular 
aperture have been corrected to total magnitude by apply- 
ing a correction of 0.4 magnitudes. Colors in BVIZYH have 
been measured running SExtractor in dual image mode, 



Grazian et al.: A critical analysis of the z ~ 7 UV LF 



3 



3 








Y-J 



Fig. 1. The Z — Y vs Y — J color-color diagram for the 
HUDF dataset limited at J < 29.2. Big points are ob- 
served galaxies, small black points are CB07 models at 
0.0 < E(B — V) < 0.3, with small red points those at 
z > 6.5. The small blue points are the result of our simu- 
lations, with synthetic objects inserted in the FITS images 
and recovered using the same procedure used for the ob- 
served catalog. The solid line is the color cut of the A crite- 
rion, while the dashed line shows the B criterion (described 
in Sect. 3). Candidate galaxies at z ~ 7 selected with cri- 
terion C are shown in magenta. Vertical arrows are upper 
limits in the Z band, while the blue circle at Z — Y > 3.6 
is the suspected SNa outlined by IQesch et al. 20101 

using isophotal magnitudes for all the galaxies, indepen- 
dently from their area. Since the FWHM of the ACS im- 
ages (0.12") is better than the WFC3 ones, we smoothed 
the ACS bands with an appropriate kernel to reproduce the 
resolution of the NIR WFC3 images. This ensures both pre- 
cise colors for extended objects and non-biased photometry 
for faint sources. 

As a comparison, we produced, only for the HUDF field, 
an equivalent BVIZYJH multiwavclength catalog, using the 
SExtractor aperture magnitudes instead of the isophotal 
magnitudes to derive galaxy colors. In the resulting Z-Y 
vs Y-J plot, indeed, galaxies are much more scattered than 
in Fig[TJ and the bulk of low-z galaxies are scattered to 
Z — Y ~ 1 when the aperture colors have been used, while 
it is limited to Z = Y ~ 0.5 if isophotal colors are adopted. 
This is due to the fact that the resulting isophotal areas are 
smaller than the circular apertures adopted, and thus are 
less noisy. 

3. Selection of z=7 candidates: color criteria 

The selection of galaxies at z ~ 7 uses the well known 
"drop-out" or "Lyman-break" technique. At 6.5 < z < 7.5, 
this feature is sampled by the large Z — Y color, as shown 
in Figure □ for the HUDF field. 



Recently, a number of slightly different color criteria 
have been used to select bona fide galaxies at z ~ 7, 
rejecting lower-z objects and cold stars entering the z- 
dropout criteria due to photometric scatter. Despite the 
differences in these criteria, however, the results obtained 
i.e. by C10a,C10blWilkins et al. 20101 IBunker et al. 20101 
IQesch et al. 20101 are comparable, at least at the bright lu- 
minosities. We aim here at applying the different color crite- 
ria to check whether the LF, which is derived from the data 
after detailed simulations, is robust against all the random 
and systematic effects that can affect each selection crite- 
rion. 

With WFC3 data (ERS and HUDF), we explo re three 
different color criteria, one of IQesch et al. 20101 one of 
IBunker et al. 20101 and one proposed by us: 

— A) The IQesch et al. 20101 criterion for the HUDF is: 

z - Y 105 > 0.8, 

z-Y W5 > 0.9 + 0.75(Yio 5 - J125), 

z-Y 105 > -1.1 + 4.0(Y 105 - J125), 

S/N(J 125 ) > 5, S/N{Y W5 ) > 5, 

S/N(V) < 2, S/N(I) < 2. 

For the ERS, not analysed in IQesch et al. 20101 an 
equivalent criterion has been adopted here, that takes 
into account the difference between the Yo9g filter in- 
stead of the Y105 one: 

z — Y398 > 1-1) 

z - Y 098 > 0.55 + 1.25(F 98 - J125), 

z-Y 098 > -0.5 + 2.0(Y 098 - J 125 ), 

S/N{J 125 ) > 5, S/N(Y 098 ) > 5, 

S/N(V) < 2, S/N(I) < 2. 

For the HUDF field, the S/N(J) > 5 criterion corre- 
sponds to J < 29.0, while for the ERS field it allows 
the selection of galaxies brighter than J = 27. 4. 

— B) The selection criterion of IBunker et al. 20101 is: 

z-Y > 1.3, 
S/N(V) < 2, 
S/N(I) < 2. 

The samples of Bunker ct al. 20101 and 
IWilkins et al. 20101 IWilkins et al. 20111 are limited 
to Y = 28.5 for the HUDF and Y = 27.2 for the ERS 
field, respectively, aiming at a more robust selection of 
z ~ 7 candidates. We adopt the same magnitude limits 
in this case. 

— C) We propose a modification of the IQesch et al. 2 010 
criterion, based on simple considerations. We adopt the 
same color criteria: 

z - Y 105 > 0.8, 

z-Y 105 > 0.9 + 0.75(Yio5 - J125), 
z-Y W5 > -i.i + 4.0(Yio5- J125) 
for the HUDF and 

z ^ Y)98 > 1-1, 

z-Y 09S > 0.55 + 1.25(^098 - Jjm), 
z-Y 098 > -0.5 + 2.0(Y 098 - J 125 ) 
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for the ERS. 

To avoid contamination of faint spurious sources, we 
add the criterion J125 — -H160 > —1.0. This additional 
criterion removes only two extremely faint galaxies with 
29.0 < J125 < 29.2 magnitudes, which are detected only 
in J but not in the H band. Indeed, the observed J — 
H color for our z-dropout candidates shows a bimodal 
distribution, with a main peak at J125 — -Hi60 ~ —0.2 
and two isolated sources at J125 — -H160 ~ —1.2. Since 
LBGs are expected to be almost flat at A ~ 2000A rest 
frame, we interpret these source with J125— i?i60 < — 1.0 
as spurious contaminants and remove them from the 
selected sample. As we have verified in the following 
(see section 5.1), this further selection does not change 
the resulting luminosity function at z ~ 7. 
For the non-detection in bands bluer than Z, we 
adopt the same criteria used in lCastellano et aL~2 010a, 
ICastellano et al. 2010bl {S/N < 2 in all BVI bands and 
S/N < 1 in at least two of them). We apply however 
only a single cut in the J magnitude (J < 29.2 for the 
HUDF and J < 27.4 for the ERS), since this band cor- 
responds to 1500A rest frame. We avoid the double 
selection S/N{J) > 5 and S/N(Y) > 5, as adopted by 
IQesch et al. 20101 since it could be less sensitive to red 
galaxies in the Y — J color at faint J magnitudes. 

FigH] shows the Z-Y vs Y-J color plot for the HUDF 
field. The solid segments show the A and C criteria, while 
the criterion B is indicated by the dashed line. 

We select 13, 6, and 20 candidates following the A, 
B, and C criteria in the HUDF field, while in the ERS 
we find 14, 6, and 22 z-dropout candidates, respectively. 
The small numbers for the candidates provided by the 
B criterion reflects the more conservative cut adopted by 
Bunker et al. 2010| to get rid of the lower-z contaminants 
with z — Y ~ 1. We have excluded from these candi- 
dates the bright transient (probably a SNa), outlined by 
IQesch et al. 20101 in the HUDF field (blue circle in FigH] at 
Z — Y > 3.6). This relatively bright transient was also iden- 
tified by iMcLure et al. 20101 and IBunker et al. 20101 while 
lYan et al. 201 II showed its possible light curve. 

Table[T]and[2]summarize the properties of the z-dropout 
candidates in the HUDF and ERS field, respectively. 

The number counts of the selected z > 7 candidates are 
different for the three different criteria. In the next para- 
graphs we will check if the simulations, by including the 
different selection effects, are able to correct for these dif- 
ferences, and if the resulting Luminosity Function depends 
on the color criteria adopted. 



4. Simulations 

While the selection criteria described above are formally de- 
signed to select a pure sample of high-z candidates, they are 
in practice applied to very faint objects, typically close to 
the limiting depth of the images. At these limits, systemat- 
ics may significantly affect their detection and the accurate 
estimate of their large color terms. To take into account all 
the systematic effects (completeness, photometric scatter) 
involved in the LF estimate, we carried out a set of detailed 
simulations, with the aim of deriving the LF. 



4.1. The Stepwise method to compute the LF 

This method is similar to that used in lBouwens et al.~ 2008, 
improved to take into account the photometric scatter, 
as described in detail in ClOa and ClOb. We first pro- 
duce a set of (Mi5oo,z) according to a flat distribution 
both in redshift and in absolute magnitudes, spanning the 
redshift range 5.5 < z < 8.5 and the absolute magni- 
tude interval —25 < M1500 < —16. We then convert the 
absolute magnitudes into observed magnitudes convolving 
the synthetic SEDs of the LBGs through the relevant fil- 
ter respo nse curves. Fo r this purpose we have used the 
models of IBruzual 20071 (hereafter CB07) with the follow- 
ing range of free parameters: metallicity: 0.02, 0.2, and 
1 Zq] age from 0.01 Gyr to the maximal age of the 
Universe at a given z; E(B-V) from to 0.3 following a 
ICalzetti et al. 20001 extinction law. Since Lyman-a emis- 
sion has an important influence on the selection function 
of Lyman break galaxies, as shown in lower redshift sam- 
ples (e.g. |Stanway et al. 2008[ |Dow-Hygehmd et al. 2007| ), 
we explicitly take into account its effect by considering a 
Gaussian distribution of Ly-a rest-frame equivalent width 
with FWHM of 30 A and centered at EW=0A. We have 
also added the intergalactic absorption using the average 
evolution as in lMadau 19951 

These synthetic galaxies are finally placed at random 
positions on the real images. We then recover the mea- 
sured magnitudes in the synthetic images using SExtractor 
with the same parameters adopted for the real images. We 
simulate all the available bands, i.e. from the B to the H 
bands for the ERS and HUDF fields. To avoid an excessive 
and unphysical crowding in the simulated images, we have 
included only 200 objects of the same flux and morphology 
each time, after masking the regions of the images where 
real objects have been detected. We repeated the simula- 
tion until a total of at least 2 x 10 4 objects were tested 
for each WFC3 field and for each morphological template 
adopted. These simulations provide the transfer function 
between the input absolute magnitudes and the observed 
colors. In addition, the simulations can be used to evaluate 
the uncertainties in the estimate of the color criteria that 
we shall exploit to detect z > 6.5 candidates, especially in 
the Z — Y and Y — J combination, which is essential for 
a clear identification of z ~ 7 galaxies. These simulations 
indeed can be used to estimate the impact of different sys- 
tematics in the study of the LF. For example, since the Y 
band filter used in the ERS field (Y 09S ) has a bluer effec- 
tive wavelength than the one in the HUDF field (Y105), it 
turns out that the selection function in redshift in these two 
fields is different. These systematics have been taken fully 
into account with these simulations. 

The simulations described above have been used to de- 
rive an estimate of the z ~ 7 LF with the Stepwise method. 
This method assumes that the rest-frame luminosity func- 
tion of galaxies can be approximated by a binned distri- 
bution, where the number density fa in each bin is a free 
parameter: this non-parametric approach allows us to con- 
strain the number density of galaxies at different magni- 
tudes without assuming an a priori functional form (i.e. a 
Schechter-like shape). 

We have assumed that the LF is made of four bins in 
the interval —22.5 < M1500 < —17.0, corresponding to the 
range sampled by our observations. We also assume that 
galaxies are uniformly distributed within the bins, with 
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-27:46:51. 


.29 


28.12 


1. 


.11 


-0.14 


-0.10 


0. 


.133 


A,C 


6,8,10,11,13,14 


HUDF- 


■2020 


03: 


:32:39. 


.725 


-27:46:21. 


.34 


28.19 


2. 


.20 


0.21 


-0.15 


0. 


.190 


A,C 


6,8,9,10,11,14 


HUDF- 


-421 


03: 


:32:37. 


.796 


-27:47:40 


.44 


28.21 


0. 


.92 


-0.01 


0.14 


0. 


.169 


C 


6,8 


HUDF- 


-2347 


03: 


:32:41. 


.824 


-27:46:11 


.26 


28.34 


1. 


.16 


0.27 


-0.17 


0. 


.180 


c 


10,11 


HUDF- 


-805 


03: 


:32:40. 


.567 


-27:46:43 


.55 


28.34 


1. 


.18 


-0.14 


-0.22 


0. 


.148 


c 


6,8,9,10,11,13,14 


HUDF- 


-1276 


03: 


:32:41. 


.596 


-27:47:04 


.47 


28.49 


1. 


.18 


-0.23 


-0.15 


0. 


.152 


c 


14 


HUDF- 


-1102 


03: 


:32:33. 


.130 


-27:46:54. 


.47 


28.73 


1. 


.57 


0.64 


-0.16 


0. 


.129 


c 




HUDF- 


-133 


03: 


:32:40. 


.334 


-27:48:02 


.61 


29.03 


1. 


.61 


0.62 


-0.15 


0. 


.109 


c 





References: [1] Bouwens et al. (2004), [2] Bouwens & Illingworth (2006), [3] Labbe et al. (2006), [4] Bouwens et al. (2008), [5] 
Oesch et al. (2009), [6] Oesch et al. (2010), [7] Fontana et al. (2010), [8] McLure et al. (2010), [9] Bunker et al. (2010), [10] Yan et 
al. (2010), [11] Finkelstein et al. (2010), [12] Castellano et al. (2010a), [13] Wilkins et al. (2011), [14] Bouwens et al. (2011). 
Note (a): both [8] and [14] indicate a nearby object 03:32:44.74 -27:46:44.92. 



Table 2. Galaxy candidates at z ~ 7 in the ERS field 



id 






RAD 


DEC 




J 


z- 


-Y 


Y-J 


J-H 




rh 


Sel 


References 


ERS- 


-5659 


03: 


:32:22.658 


-27:43:00. 


.64 


25.54 


1. 


.34 


0.25 


-0.01 


0. 


202 


A,B,C 


1,2,3,5 


ERS- 


-3881 


03: 


:32:25.284 


-27:43:24. 


.20 


25.98 


2. 


.07 


0.11 


0.14 


0. 


.254 


A,B,C 


4,5 


ERS- 


-7034 


03: 


:32:24.095 


-27:42:13. 


91 


26.08 


1. 


.22 


0.38 


-0.17 


0. 


.275 


C 


4,5,6 


ERS- 


-285 


03: 


:32:10.407 


-27:45:40. 


.80 


26.43 


1. 


.13 


-0.01 


0.17 


0. 


.225 


C 




ERS- 


-1819 


03: 


:32:06.826 


-27:44:22. 


18 


26.61 


1. 


.26 


-0.07 


-0.30 


0. 


135 


A,C 


6 


ERS- 


-2200 


03: 


:32:22.934 


-27:44:09. 


92 


26.72 


1. 


.75 


0.10 


0.01 


0. 


.128 


A,B,C 


5 


ERS- 


-3719 


03: 


:32:15.400 


-27:43:28. 


61 


26.98 


1. 


.31 


-0.08 


-0.33 


0. 


.120 


A,B,C 




ERS- 


-4682 


03: 


:32:23.155 


-27:42:04. 


69 


27.00 


1. 


.88 


0.09 


-0.61 


0. 


.178 


A,C 


5 


ERS- 


-10069 


03: 


:32:22.518 


-27:41:17. 


34 


27.01 


1. 


.22 


-0.01 


-0.10 


0. 


.129 


A,C 


6 


ERS- 


-1676 


03: 


32:08.359 


-27:44:27. 


40 


27.03 


1. 


.95 


0.29 


0.01 


0. 


159 


A,C 




ERS- 


-457 


03: 


:32:10.035 


-27:45:24. 


53 


27.09 


3. 


.70 


0.02 


-0.47 


0. 


.146 


A,C 




ERS- 


-8656 


03: 


:32:36.313 


-27:40:14. 


99 


27.11 


2. 


.34 


0.12 


0.07 


0. 


.197 


A,C 




ERS- 


-3875 


03: 


:32:09.853 


-27:43:24. 


01 


27.13 


3. 


.97 


-0.29 


0.21 


0. 


.185 


A,B,C 




ERS- 


-6532 


03: 


:32:31.227 


-27:42:25. 


17 


27.14 


1. 


.86 


-0.25 


-0.88 


0. 


136 


A,B,C 




ERS- 


-3642 


03: 


:32:41.793 


-27:43:30. 


05 


27.16 


1. 


.59 


0.46 


-0.12 


0. 


357 


A,C 




ERS- 


-4752 


03: 


:32:39.829 


-27:42:40. 


90 


27.19 


1. 


.69 


0.67 


-0.38 


0. 


.288 


C 




ERS- 


-7187 


03: 


:32:33.860 


-27:42:09. 


73 


27.30 


3. 


.19 


0.34 


-0.11 


0. 


.140 


C 




ERS- 


-148 


03: 


:32:07.287 


-27:45:54. 


29 


27.31 


2. 


.90 


-0.18 


-0.63 


0. 


130 


A,C 




ERS- 


-4639 


03: 


:32:13.405 


-27:42:30. 


89 


27.33 


2. 


.50 


0.77 


-0.02 


0. 


.168 


C 




ERS- 


-5945 


03: 


:32:16.767 


-27:43:07. 


16 


27.35 


3. 


.02 


0.25 


-0.22 


0. 


.145 


C 




ERS- 


-5383 


03: 


:32:29.458 


-27:42:54. 


27 


27.36 


2. 


.09 


0.54 


-0.66 


0. 


136 


c 




ERS- 


-1965 


03: 


:32:14.988 


-27:44:17. 


54 


27.37 


2. 


.82 


0.04 


-0.67 


0. 


.148 


c 





References: [1] Bouwens et al. (2010b), [2] Hickey et al. (2010), [3] Castellano et al. (2010a), [4] Wilkins et al. (2010), [5] Wilkins 

et al. (2011), [6] Bouwens et al. (2011). 



number densities 4>i to be determined. With these simulated 
galaxies, we compute the distribution of observed magni- 
tudes originated from each bin, scaled to the observed area 
in the WFC3 fields (subtracting the area covered by de- 
tected objects). We then find the combination of fa that 
best reproduces the magnitude distribution of our observed 
objects with a simple x 2 minimization. Since the photo- 
metric scatter can move galaxies from a magnitude bin to 
the near ones, the \ 2 minimization ensures that this ef- 
fect is taken into account. The resulting uncertainties can 
be higher than those derived assuming a simple linear re- 



lation between the apparent and absolute magnitude, as 
in iBouwcns et al. 2008] These differences are discussed in 
detail in [Castella no et al. 2010bl 

Finally, the simulations are used to estimate the sys- 
tematic effects acting when we use colors at shorter wave- 
lengths, i.e. in the BVI bands, to reject possible interlop- 
ers. Because of the large IGM and internal HI absorption, 
the expected flux in these bands for z > 6.5 galaxies is far 
below the detection threshold, or even null. For this rea- 
son a stringent limit on the measured flux in these bands 
is adopted to remove lower redshift interlopers. However, 
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the S/N estimated by SExtractor may be a poor represen- 
tation of the actual photometric scatter at low fluxes, due 
to a combination of factors, such as uncertainties in the 
estimate of the local background, underestimates of the 
true r.m.s., or chance superposition of faint blue galaxies 
along the line of sight. To account for these effects, we have 
measured the resulting signal-to-noise (S/N) ratio in the 
BVI images for each simulated object inserted in the J 
one, which should be zero on average. It turns out that 
the actual distribution of the S/N ratios is wider than the 
one obtained with SExtractor, which is computed scaling 
the input weight image. We thus estimate the "effective" 
r.m.s. cf s/N: the r.m.s. of the signal-to-noise distribu- 
tion in each of the BVI images, which is about 1.5 for the 
B and V bands and 1.8 for the I band. Even taking into 
account this wider distribution, we find that the tails of 
the S/N distribution (S/N > 2ag/ N ) contain more objects 
than in the case of a pure Gaussian distribution. As men- 
tioned above, we will use the estimated <js/n in all BVI 
bands, requiring that high-z candidates have flux < 2<js/n 
in all BVI bands and < lcrg/N in at least two of them. With 
our simulation, we estimate that the fraction of true high-z 
galaxies lost because of this strict criterion is about ~ 12% 
for the HUDF field. For the ERS, we find an effective S/N 
ratio of 1.8 for all the BVI bands and a rejection of true 
high-z galaxies of ~ 24%. Such an aggressive requirement 
is needed in order to clean our sample by low-z interlopers 
that are contaminating our z-dropout sample. This tech- 
nique has been extensively discussed in ClOa and CI Ob. In 
particular, if we neglect the requirement that S/N must be 
< la s/n in at least two bands out of the BVI ones, we find 
32 and 97 galaxy candidates for the HUDF and ERS fields, 
respectively, to be compared with our 20 and 22 master can- 
didates. We thus prefer to reduce significantly the numbers 
of possible contamination by interlopers, at the expenses of 
being less complete by 24%, especially for the ERS field. 
This fraction is significantly higher than that found in the 
HUDF and it is probably due to the fact that the depth of 
the ACS GOODS images is not matched to the deep lim- 
iting magnitude one can reach with WFC3 in the ERS. In 
this sense, these simulations are essential to estimate the 
LF at z ~ 7. 



4.2. The morphological profiles used during simulations 

We expect that the output of the simulations depends criti- 
cally on the assumed morphology. As realistic morphologies 
for the simulations related to the WFC3 fields, we check 
different profiles. In particular, we explored the following 
morphologies spanning from the most compact to the most 
extended objects available in the WFC3 fields, exploring 
also plausible extrapolations from lower-z studies: 

— 1) We select four point- like objects in the HUDF field 
based on their FWHM, avoiding bright saturated ob- 
jects. The spectroscopic classification is not available 
for these objects, and their colors are consistent with 
those of stars. Since these are compact, we use them as 
an input profile for very compact galaxies at z = 7, as it 
is expected that faint galaxies at high-z are very small 
(j Windhorst et al 2002]) . 

- 2) The transient object at Z - Y > 3.6 in FigQ] (prob- 
ably a SNa) is used as an empirical profile for our sim- 
ulations. Since this object was not detected in previ- 



ous NICMOS images of HUDF (IQesch et al. 20101) . it 
is plausible that the host galaxy in the background is 
not able to change the SNa profile, and the result is a 
very compact object. We used this template to check 
whether it gives different results from the four stars de- 
scribed above. 

— 3) We select the three brightest and most compact z ~ 7 
candidates in the HUDF (J < 27.4) as input for our 
simulations, assuming that the physical dimension of 
high-z galaxies does not vary dramatically with their 
magnitudes. We avoid in this group the object HUDF- 
1095 (G2.1408 of ClOa) since it is the most extended 
z ~ 7 candidate in our sample and we consider it in a 
separate group. 

— 4) We study the effect of very extended objects or of 
merging at high-z using as input in our simulations the 
object HUDF-1095 (G2.1408 of ClOa). It is very bright 
with two knots, possibly the result of a merging of two 
galaxies. Alternatively, its clumpy morphology in the 
rest frame UV could be due to extended star forming 
knots in a nascent galaxy. First results based on spectro- 
scopic confirmation with FORS2 seem to indicate that 
it is a z = 6.972 LBG (|Fontana et al. 2010|) . confirming 
the photometric redshifts from ACS and WFC3 based 
SED. Though not representative of the whole z = 7 
population, this galaxy is an LBG at z ~ 7 and it is in- 
teresting to study the resulting LF assuming this mor- 
phology as input for our simulations. 

— 5) There is the possibility that fainter galaxies are more 
compact than brighter ones, as expected in the hierar- 
chical formation models. To check this hypothesis, we 
adopt as input profile the one obtained by stacking the 
15 fainter galaxy candidates at z=7 (27.5 < J < 29.2) 
in the HUDF field. 

— 6) We used as input for our simulations the pro- 
files of spectroscopically confirmed galaxies at z ~ 
4 in the ACS I-band of the HUDF, scaled as 
(1 + z) _1 in their physical sizes, as adopted in 
IBouwens et al. 20081 and IBouwens et al. 2010al Since 
there is a regular size evolution with redshift, 
as found by |Ferguson et al. 2004] and confirmed by 



IBouwens et al. 20041 and IQesch et al. 2010bl it is useful 
to explore this case, in order to compare our results with 
earlier LFs that have appeared in the literature. These 
templates are also useful to avoid using self-selected z=7 
candidates as input for our simulations. 

FigJU shows all the model galaxies used for the HUDF 
and ERS simulations. 

In addition of adopting single objects for the simula- 
tions, we explore also different distributions of half light 
radii (rh) at z ~ 7: 



— 7) Following Oesc h et al. 20071 we explore a uniform 
distribution in the half light radius of the synthetic 
galaxies, assuming an exponential declining profile. The 
two extremes of the uniform distribution are derived 
from the values used by IQesch et al. 20071 for their sim- 
ulations to recover the z=5 LF from HUDF data. They 
adopt a uniform distribution with half light radius from 
0.05 to 0.5 arcsec. Converting these values to physical 
units at z=5, and applying a (1 + z) -1 evolution until 
z=7, and reconverting in angular dimensions, we derive 
the limits to be 0.045 to 0.45 arcsec. Based on FigOJ it 
is clear that this distribution can be ruled out by the 
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Fig. 2. The different morphologies used for the simula- 
tions of WFC3 fields. From top left, the first four cutouts 
represent non-saturated stars in the HUDF survey (mor- 
phological template 1). The SNa candidate outlined by 
IQesch et al. 20101 (template 2) is then shown, followed by 
the three brightest z ~ 7 candidates of the HUDF sur- 
vey (template 3). The most extended galaxy is the bright- 
est candidate in the H UDF (HUDF-1095 or G2.1408 of 
ICastellano et al. 2010al template 4) . The second thumbnail 
on the third raw shows the stack of the faint candidates se- 
lected in the HUDF (template 5). All these thumbnails have 
been extracted from the J125 band of HUDF. The remain- 
ing objects are spectroscopically confirmed galaxies z ~ 4, 
extracted from the ACS-I band of the HUDF and shifted at 
z=7 as described in the main text (template 6), and two ar- 
tificial exponential disk galaxies with half light radius=0.3 
arcsec (templates 7 and 8) and with half light radius=0.1 
arcsec (templates 9 and 10). 



present observations. However, we decided to study this 
uniform distribution as a "thought experiment" on how 
hypothetical size distributions might change the derived 
completeness at faint magnitudes, but we do not use it 
to estimate the z ~ 7 luminosity function. 



— 8) Another possible distribution for the half light ra- 
dius of the simulated galaxies is the log-normal distri- 
bution with mean and dispersion adopted following the 
relation of [Fterguson et al. 2004 for photometrically se- 
lected LBGs at z=4, applying a (1 + evolution to 
the half light radii in physical units. We derive for the 
z=7 population a mean half light radius of 0.258 arc- 
sec and a = 0.0358. This approach has been adopted 
also in IQesch et al. 20091 and IQesch et al. 20101 In the 
following it is indicated also as Ml. 

— 9) The typical half light radii of z=7 candidates are of 
the order of 0.1-0.15 arcsec. As in IQesch et al. 2 010b, 
we assume a log-normal distribution of the half light 
radius of the simulated galaxies with mean=0.153 and 
sigma=0.0574 arcsec. This distribution, however, can- 
not explain the presence of the HUDF-1095 galaxy, with 
rh > 0.25arcsec. In the following it is indicated also as 
M2. 

— 10) A different distribution for the half light radius 
of the simulated galaxies at z ~ 7 has been derived 
starting from the observed half light radii of spectro- 
scopically confirmed LBGs at z = 4 in GOODS from 
IVanzella et al. 20091 We select only galaxies with 24 < 
I < 25 magnitudes to avoid incompleteness for extended 
LBGs. Applying a (1 + z) _1 evolution to the half light 
radii at z ~ 4 in physical units, the log-normal distri- 
bution has a mean=0.177 and sigma=0.0051 arcsec at 
z ~ 7 and in apparent dimensions. In the following it is 
indicated also as M3. 

From FigH] it is clear that the observed distribution of 
half light radii for the z=7 candidates peaks at ~ 0.1 arcsec- 
onds and extends to 0.2-0.25 arcsec. The small area covered 
by the two WFC3 surveys considered here is at present too 
small to sample the extended tail of the half light radii 
distribution, and it is not possible to infer a robust size dis- 
tribution. Clearly, the uniform distribution (criterion 7) is 
not representative of the observed z = 7 galaxy population 
and will not be considered in the following to estimate a 
reliable z ~ 7 luminosity function. 

The three log-normal distributions (templates 8,9, and 
10 shown in FigJ3]) cannot be confirmed or ruled out at the 
present stage for a number of reasons: 

— a) the z ~ 7 candidate sample at magnitude J < 27.8 
in the HUDF is limited to 6 galaxies. We have chosen 
this magnitude limit since it is known from simulations 
that in the HUDF we are complete at the ~ 80% level 
for extended galaxies (see FigJH blue and cyan lines). 
The corresponding limits for shallower surveys (ERS or 
the parallel HUDF observations HUDF-P1 and HUDF- 
P2) are 1-2 magnitude brighter (J < 26.0 - 26.8) than 
in the HUDF. Adopting these limits, only two addi- 
tional brighter candidates have been found on the ERS 
survey, both having rh > 0.2arcsec. Since the number 
statistic here is very low due to the limited area covered 
by ultradeep HST observations, at the present stage we 
cannot exclude the presence of more extended galaxies 
at J < 27.8 magnitude. 

— b) the extrapolation to fainter magnitudes of the size 
distribution for z ~ 7 galaxies is not known at present. 
It is possible that a population of extended sources (sim- 
ilar to HUDF-1095) exists at high-z and at fainter mag- 
nitudes. The present HUDF observations cannot rule 
out such hypothesis. 
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Fig. 3. The black histogram shows the observed distri- 
bution of the half light radii (in arcsec) of the z ~ 7 
galaxy candidates brighter than J = 27.8 in the HUDF. 
At this magnitude limit the completeness of the J band 
image is higher than 80% even for the most extended 
objects. The galaxy with half light radius ~ 0.25 arcsec 
is the HUDF-1095 one (template 4). The cyan solid line 
shows the uniform distribution in half light radius (tem- 
plate 7), while the dark pink and green solid lines in- 
dicate the two log-normal distributions in half light ra- 
dius adopted by Ferguson et al. 2004] (template 8, Ml) 
and bv lOesch et al. 2010bl (template 9. M2), respectively. 
The purple line indicates the distribution derived from 
IVanzella et al. 20091 (template 10, M3). All sizes are cor- 
rected for PSF. The dashed line shows the completeness of 
galaxy detection for different half-light radii derived by the 
performed simulations. This curve has been multiplied by 
a factor of 2 only for clarity. 



c) the spectroscopic confirmation of z ~ 7 galaxy candi- 
dates with detailed HST morphology and half light ra- 
dius determination in the HUDF and ERS fields is still 
lacking at this stage and some of them could be lower- 
z contaminants. The brightest galaxy in the HUDF 
(HUDF-1095, which is clumpy or in a merging state) 
turned out to be spectroscopically confirmed at z=6.972 
(|Fontana et al. 2010|) . Since it is the only galaxy with 
detailed HST morphology in the NUV rest frame (tem- 
plate 4) to be confirmed at z ~ 7 in the HUDF field, 
and it is relatively bright (J = 26.07) and extended 
(rh = 0.256), we cannot exclude that a similar popula- 
tion of extended sources exists at high-z and at fainter 
magnitudes . 

d) in IBouwens et al. 20111 and also in 
Wilk ins et al. 20111 a z ~ 7 galaxy candidate with 
an half-light radius of 0.46 arcsec has been found 
(UDF092z-01191133). Moreover, other z ~ 7 candi- 
dates found in the ERS, HUDF-P1, and HUDF-P2 
fields show extended half light radii (0.3 arcsec or 



more). If they turn out to be spectroscopically con- 
firmed at z ~ 7, then distributions with extended half 
light radii, similar to template 8 (Ml), would become 
more reliable. 

Based on these considerations, it is not possible at the 
present stage to exclude any of the adopted log-normal half 
light radii distributions (Ml, M2 or M3). 

The simulations carried out with the different morpho- 
logical templates described above have been used to esti- 
mate the completeness in the detection procedure. This is 
shown in Figure El where we plot the fraction of detected 
versus input objects as a function of the J125 input magni- 
tude for the HUDF field. 

The completeness for the HUDF is basically equal to 
100% down to J ~ 27 for all the adopted morphologies, and 
fades in different ways according to the different half light 
radii of the templates adopted in the simulation. In partic- 
ular, more extended galaxies typically have a lower com- 
pleteness at J > 27, and the fraction of faint z—7 galaxies 
recovered depends dramatically on the morphological tem- 
plates used during the simulations. Since the completeness 
is used to estimate the LF, we expect the faint end of the 
z=7 LF to be strongly dependent on the morphological tem- 
plate adopted, as we will show in the next section. 

5. Results 

The derived LF is sensitive to the details of the simulations 
adopted. Here we are interested in exploring mainly two 
issues: the role of the color selection criteria on the esti- 
mate of the LF and the role of galaxy morphology on the 
measurement of a, the steepness of the faint end of the LF. 

5.1. Results for different color selections 

We first explore three different color criteria for the se- 
lection of z ~ 7 galaxies. The LF found with these three 
selection criteria (A, B, and C), described in Section 3, 
for the HUDF and ERS fields is shown in FigE] The red 
filled circles show the LF derived with the color criterion C, 
compared to the one derived with the B (small filled black 
triangles) and A criteria (blue squares), for the HUDF data 
alone. We plot for comparison also the Hawk-I results (ma- 
genta open circle) of ClOa. The WFC3 LF has been derived 
assuming as morphological inputs for our simulations real 
bright (J < 27.5) candidates at z ~ 7 found in the HUDF 
(template 3). Adopting as morphological templates spec- 
troscopically confirmed galaxies at z ~ 4 found in the ACS 
i-band of the HUDF and scaled to z ~ 7 (template 6), or 
the stack of the fainter (J > 27.5) candidates (template 5), 
we obtained very similar results. We have also verified that 
the additional limit J125 — -H160 > —1.0 adopted in crite- 
rion C does not imply significant changes in the LF in the 
faintest bin {M uv ~ -18.5). 

It is clear that the three color selection criteria analysed 
here give comparable results, well within the uncertainties, 
despite the size of the three samples changes by a factor of 
^3. This is due to the fact that all the three criteria are well 
designed to isolate the few z ~ 7 candidates from the great 
number of interlopers, but each method makes different 
assumptions, resulting in a variable number of z-dropout 
galaxies. The resulting LFs, however, are similar, since the 
simulations reguired to derive them are able to take into 
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Fig. 4. The fraction of detected galaxies versus the J125 
input magnitude for the simulations in the HUDF field. 
Different colors indicate the various morphologies adopted: 
the dashed magenta lines indicate non-saturated stars on 
the HUDF survey (morphological template 1). The solid 
magenta line shows the completeness obtained using the 
SNa candidate outlined by IQesch et al. 20101 (template 2) 
as morphological template, while the red solid curves show 
the three brightest z ~ 7 candidates of the HUDF survey 
(template 3). The blue solid line is the completeness de- 
rived by adopting as template the brightest candidate in the 
HUDF (HUDF-1095 or G2.1408 of lCastellano et al. 2010al 
template 4), and the black solid line represents the com- 
pleteness obtained adopting as input template the stack of 
the faint candidates selected in the HUDF (template 5). 
The green dashed lines indicate the spectroscopically con- 
firmed galaxies z ~ 4, shifted at z = 7 (template 6). The 
cyan solid line shows the completeness obtained assuming a 
uniform distribution in half light radius (template 7) , while 
the dark pink and green solid lines indicates the results of 
two log-normal distributions in half light radius adopted by 
|Ferguson et al. 2004| IQesch et al. 20071 lOesch et al. 20101 
(template 8, Ml) and by lOesch et al. 2010bl (template 9, 
M2), respectively. The purple line ind icates the complete- 
ness derived from IVanzella et al. 20091 (template 10, M3). 



account properly the random and systematics effects (in- 
completeness and contamination by photometric scatter), 
especially at the faintest magnitudes. In the faintest magni- 
tude bin the sizes of the error bars depend on the criteria 
adopted, with the larger errors for the most strict criterion. 
The B criterion selects indeed a smaller number of candi- 
dates, with a resulting higher Poisson noise than the one 
derived adopting the A or C criteria. 

The only difference between the A and C criteria is 
the magnitude cut adopted: we s elect all galaxies down to 
J = 29.2, while lOesch et al. 2010l recover only galaxies with 
S/N greater than 5 in the Y and J bands. These color cri- 
teria are similar, and no systematic differences are found in 



Fig. 5. The LF at z=7, obtained using as morphological 
templates real bright (J < 27.4) candidates at z ~ 7 found 
in the HUDF (template 3). The magenta circle represents 
the LF found in the Hawk-I from ClOa, while the other 
symbols indicate the LF in the HUDF and ERS derived 
adopting different color criteria: our color criteria (red filled 
circles, criterion C), IBunker et al. 20101 (small filled black 
triangles, criterion B) and lOesch et al. 20101 (blue squares, 
criterion A). The blue points are shifted by -0.1 magni- 
tudes and the black points by +0.1 for clarity. The blue 
dashed and green dot-dashed lines show the LF at z=6 of 
IMcLure et al. 200"9l and iBouwens et al. 20071 respectively. 



the LF estimate. For this reason, we decide to adopt the C 
color criterion in the following sections. 



5.2. Results for different morphologies 

We then explore the issue of adopting different morpholog- 
ical templates in the simulations. For each morphological 
type shown in FigJ5] and discussed in Section 4, we have 
simulated 2 x 10 4 synthetic galaxies as described in Sect. 
4 and used them to compute the LF at z ~ 7. In Fig(5]wc 
show the LF for each template considered: at M1500 ~ — 20, 
the different morphologies adopted have a small impact on 
the number density <E>, while at M1500 ^ — 18 the LF is ex- 
tremely dependent on the assumption of the morphological 
templates. In Fig|6]the open points show the result of tem- 
plates 1-6, based on single object morphologies, with the 
relevant uncertainties. 

In Fig [7| we considered different distributions for the 
simulated half light radii, namely the uniform one (tem- 
plate 7, cyan points), and the three log-normal distribu- 
tions (template 8 using dark pink points, template 9 using 
dark green points, and template 10 using purple points). For 
comparison, we plot the LF derived by[Oesch ct al. 20101 as 
filled black squares: these are in agreement with the result 
based on the simulations carried out with the distributions 
of half light radii which extend to 0.2-0.4 arcsec. The large 



10 



Grazian et al.: A critical analysis of the z ~ 7 UV LF 




Fig. 6. The LF at z=7, obtained using as morphologi- 
cal types galaxies drawn from the templates 1-6. The color 
code is the same of Fig. 4. The blue and green bins are 
shifted by -0.2 and -0.1 magnitudes, while the red and black 
points are moved by +0.1 and +0.2, respectively, for clar- 
ity. The blue dotted lines indicates the range of variability 
for the best fit of a spanning the range —2.15 < a < —1.38. 
The blue dashed and green dot-dashed lines show the LF 
at z=6 of lMcLure et al. 20091 and IBouwens et al. 20071 re- 
spectively. The fil led black squares are the LF presented in 
IQesch et al."20T0l 



variance in the galaxy number density at the faint end of 
the z ~ 7 LF depends on the different completeness cor- 
rections applied to the same dataset and thus reflects the 
uncertainties due to the poorly known half light radius dis- 
tribution. 

Fixing the Schechter parameters <J>* and M* to 6.9 x 
10~ 4 and -20.10 (the best fit values of lOuchi et al. 2009|) . 
respectively and varying a, we derive a large range of vari- 
ation for the best fit of the LF steepness, namely —2.15 < 
a < —1.38 at 68% c.L, represented by the two blue dotted 
lines in the plot. 

The steepness of the faint end of the LF, a, depends crit- 
ically on the half light radii of the synthetic galaxies used 
to carry out the simulations. In FigJS]we show the half light 
radii of the different morphological templates against the 
best fit a of the LF. For large value of the half light radius, 
namely > 0.2 arcsec, the best fit a is between -1.8 and 
-1.9, as found by IQesch et al. 20101 IBouwens et al. 2010a! 
and IBouwens et al. 20111 at z=7 and beyond. For smaller 
values, ~ 0.1 arcsec, we find a in the range -1.4 to -1.7, 
in marginal agreement with IMcLure et al. 20l"0l that have 
used point sources as input for their simulations. The un- 
certainties on this parameter are still large, of the order of 
~ 0.2 — 0.4. We find however a significant relation between 
the half light radius and a, indicating that simple conclu- 
sions based on the LF at z ~ 7 depend dramatically on the 
half light radius distribution assumed for the simulations. 



Fig. 7. The LF at z=7, obtained using different distribu- 
tions for the simulated half light radii, namely the uniform 
one (template 7, cyan points), and the three log-normal 
distributions (template 8 using dark pink points, template 
9 using dark green points, and template 10 using purple 
points). The dark green bins are shifted by +0.1 magni- 
tudes, while the dark pink points are moved by -0.1 magni- 
tudes, respectively, for clarity. The blue dotted lines indi- 
cate the range of variability for the best fit of a spanning the 
range —2.15 < a < —1.38. The blu e dashed and green dot- 
dashed lines show the LF at z=6 of IMcLure et al. 20091 and 
IBouwens et al. 20071 respectively. The filled black squares 
are the LF presented in [Oesch ct a l. 20101 

5.3. Fitting the LF at z ~ 7 

We want to derive here the LF at z ~ 7 making reliable and 
robust assumptions based on the considerations outlined 
above. We adopt the color criterion C, with a magnitude cut 
at J=27.4 and J=29.2 for the ERS and the HUDF, respec- 
tively. As shown above, the simulations are able to correct 
the incompleteness of this color criterion in a reliable way. 
The morphology templates used in the simulations, how- 
ever, are critical, since the simulations cannot correct for a 
biased assumption on the profiles of the candidates. Since 
it is not possible, with the small number statistics achieved 
by the present WFC3 data, to derive precisely the distri- 
bution of the half light radii of the z ~ 7 LBGs, especially 
at the faint end and at high values of rh, we carry out the 
fit for the two most extreme, but reliable, variants of the 
LFs derived in the previous paragraph. We have adopted 
the simulations carried out with the stacking of the faint 
candidates selected in the J band of HUDF (template 5) as 
an extreme case (all the z=7 galaxies are compact), while 
the HUDF-1095 galaxy (template 4) has been chosen as the 
opposite case (all high-z LBGs are very extended). We have 
verified that all the other options (templates 1-3 and 6-10) 
give results that are in between these two extremes. 

We have combined the results of the ERS and HUDF 
fields in a single LF, to reduce the statistical uncertain- 
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Fig. 8. The dependence of the steepness of the LF, 
a, against the half light radii of the simulated galaxies. 
Different colors indicate various morphology adopted, us- 
ing the same color code of Fig. 4. 



ties and to have an extended range in the absolute mag- 
nitude Mi 5oo- To cover the bright p art of the LF, we also 
use the results of IQuchi et al. 2 009 obtained from a shal- 
low and wide Subaru survey and those of ClOa and ClOb 
derived from an intermediate depth survey with Hawk-I. 
In Fig|9]we show our results (red and green filled sq uares) 
complemented with the results of IQuchi et al. 20091 (filled 
black circles) , Bouwen s et al. 2010al (filled cyan circles) and 
of ClOa and ClOb (empty black circles). 

All the observed points shown in Fig(3] have been 
used to fit the LF. We scan the parameter space in 3D 
($*, M*, and a) adopting a Schechter parameterization 



(Schechter(1976) ) and taking into account the asymmet- 
ric errors and the presence of upper limits. The best fit is 
shown (long dashed magenta lines), and it is compared to 
the early Hawk-I result of ClOa (black solid line) and the 
LFs at z=6 of lBouwens et al. 20071 and IMcLure et al. 20091 
Using our best fit, we provide a robust estimate of the LF at 
z = 7, excluding at >99.73% c.l. (> 3<r) the non-evolution 
from z=6 to z=7. Indeed, we are able to show that the 
normalization at A/1500 ~ —19 is significantly lower than 
the LF at z ~ 6. Moreover, combining data with a large 
range in absolute magnitudes we are able to reduce the de- 
generacies in the three parameters of the LF. In FigfTUl we 
provide the confidence regions at 68% c.l. for $*, M*, and 
a at z=7, exploring also the constraints given on the SFR 
density (SFRD). 

Table [3] summarizes the best fit parameters for the 
Schechter LFs, with their 68% c.l. uncertainties. The lu- 
minosity density puv has been computed integrating the 
fitted LFs down to an absolute magnitude (at 1500 A rest 
frame) of Mjjv = — 10 and to Mjjv = —19. We have con- 
verted these values in a SFRD following the standard for- 
mula by Madau et al. (1998) and applying the extinction 
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Fig. 9. The LF at z=7, obtained using as morphological 
templates the stacking of the faint candidates at z ~ 7 
found in the J band image of the HUDF (template 5) and 
the C color criterion. The LF is computed using ERS and 
HUDF fields together (red filled squares, G10). The green 
filled squares show the LF obtained with the same color 
criterion C but using the HUDF-1095 candidate as input 
for the simulations (template 4). These data have been 
complemented with th e results of IQuchi et al. 20091 (O09), 
IBouwens et al. 2010bl (B10) and ClOb. The lower dashed 
magenta line is the best fit to all the observed points related 
to the template 5, the upper dashed magenta line is the best 
fit related to the template 4, while the black solid line shows 
the Hawk-I result of ClOa. The blue dashed and green dot- 
dashed lines show the LF at z=6 of IMcLure et al. 2009 
(M09) and IBouwens et al. 2007 (B07), respectively. 



correction of Meurer et al. (1999) of A1600 = 0.43 (consid- 
ering an average UV slope j3 = —2.0 for the SED of the z=7 
LBG candidates). Our estimates of the luminosity density 
and SFRD are consistent with other values in the literature, 
reducing significantly their uncertainties. 



6. Discussion 

We discuss the implications on reionization derived from 
our best fit of the LF and the uncertainties on the LF pa- 
rameters implied by the assumption made during complete- 
ness simulations. 



6.1. Reionization 



Following IBolton fc Haehnelt 20071 and IQuchi et al.~2 009. 
the contribution of z=7 LBGs to the hydrogen ionizing pho- 
ton can be set as 



Nio^Mpc- 3 ) = 10 49 ' 7 ( 



PUV w 3.0 / e ,j C 

6-10 25n a s A 0.1 ' 



(1) 
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Fig. 10. The confidence regions at 68% c.l. (1 a) for the LF parameters and the SFRD. The red points show the 
projection of the best fit. The plot on the left shows the confidence regions based on the simulations with compact 
morphology (template 5), while the right plot is related to extended morphology (template 4). 



Table 3. Schechter LF fit Results 



Compact Morphology 


Best Fit 


Lower (1 a ) 


Upper(lcr) 


Log(®*) 


-3.13 


-3.66 


-2.80 


M* 


-20.14 


-20.68 


-19.74 


a 


-1.65 


-2.17 


-0.94 


puv(-W) 


9.2E+25 


3.8E+25 


7.9E+26 


SFRD(-IO) 


0.017 


0.007 


0.150 


puv(-W) 


2.7E+25 


1.9E+25 


3.5E+25 


SFRD(-19) 


0.005 


0.004 


0.007 


Extended Morphology 


Best Fit 


Lower (1 a ) 


Upper(lcr) 


Log(®*) 


-3.22 


-4.10 


-2.86 


M* 


-20.30 


-21.12 


-19.94 


a 


-2.02 


-2.61 


-1.51 


puv(-W) 


4.5E+26 


8.8E+25 


4.4E+28 


SFRD(-W) 


0.085 


0.017 


8.365 


Puv(-W) 


3.2E+25 


2.3E+25 


3.9E+25 


SFRD(-19) 


0.006 


0.004 


0.007 



puv(-W) is the luminosity density obtained integrated the LF 



down to an absolute magnitude of -10 based on a simple 
extrapolation, while puv{— 19) is the one derived integrating 
the best fit LF down to an absolute magnitude M1500 of -19. 
The SFRD(-W) and SFRD (-19) represent the SFR density 

limited to -10 and -19 in Misoo- The unity is in 
erg/s/Hz/Mpc 3 for puv and in M e /yr/Mpc 3 for SFRD. The 
upper part of the table summarizes the LF parameters 
resulting from simulations based on template 5 (compact 
morphology), while the lower part is related to template 4 
(extended morphology). 



where as is the spectral index of ionizing emission and 
f esc is the escape fraction of ionizing photons. The rate 
of ionizing photons needed to balance the recombination 
process of hydrogen in the IGM and hence to keep the 
Universe reionized is N ion > 10 47A Chii(1 + z) 3 , that trans- 
lates into a constraint to the dumpiness of the IGM to be 
Chii < (puv fesc)/ '(10 22 as(l + z) 3 ). We fix these parame- 
ters to z = 7, p = —2, and as — 3 in the following. In FigfTTI 
we show in the dashed areas the Chii allowed by a given 
f esc to keep the IGM ionized at z ~ 7: the red region is for a 
compact morphology (template 5), while the blue one refers 



to an extended morphology (template 4). The relevant UV 
emissivity of LBGs puv is computed by integra ting the 
LF down to Mi 5 o = — 10 (jBouwens et al. 2011|) . assum- 
ing that the steepness of the faint end of the LF remains 
constant; we have no information on the number density of 
z ~ 7 galaxies at these faint magnitudes with the present 
data. Given our confidence intervals for the LF parame- 
ters, and hence the limit to puv(— 10), we derive a limit 
Chii < 6.5 (at 68% c.l.) in order to have the Universe ion- 
ized at z = 7, assuming a maximum escape fraction of 1.0, 
and the simulations based on template 5 (compact mor- 
phology), as shown in FigfTTI If we use the LF based on 
extended morphology (template 4), the 68% c.l. limit to 
the dumpiness of the ionized hydrogen is Chii < 29.6. In 
this reasoning we must consider Chii > 1, since galaxies at 
z ~ 7 are formed in biased density regions of the Universe 
and the IGM is plausibly not homogeneous (Chii = 1) at 
these redshifts. Limiting the dumpiness to Chii > 1, we 
thus have f esc > 0.17 in order to have the Universe reion- 
ized by z=7 assuming compact morphology, while the limit 
is fesc > 0.034 for extended morphology. 

Bolton & Haehnelt (2007) have inferred Chii < 3 at 
z ~ 6 from the Lyman-a forest photoionization state: since 
Chii is expected to monotonically decrease towards high-z 
in a hierarchical Universe, an escape fraction of f esc ~ 0.5 
is enough for the reionization of the IGM (but not for the 
more clumpy galaxy formation regions) considering tem- 
plate 5 (compact morphology), while a value of f eS c ~ 0.1 is 
sufficient according to template 4 (extended morphology). 

We neglect in our computation the contribution of 
AGNs, since their LF at z=7 is still unknown and the 
upper limits currently available indicate that the AGN 
will add only 5-8% to the luminosity density of galaxies 
(|Ouchi et al. 2009p . 

The reionization rate Ni on however depends largely on 
the contribution of faint galaxies (—19 < Afi 50 o < —10), 
since puv{^i5oa < —10) could be 10-15 times larger than 
Puv{M\5oo < —19) if the LF is steeper than a < —2. This 
is eventually still possible, if a large number of faint and 
extended galaxies exist, similar to the object HUDF-1095 
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Fig. 11. The constraints on the maximum Chii and min- 
imum f esc parameters in order to have the Universe reion- 
ized by LBGs at z=7 given our LF. Given a value for f eS c, 
the Universe is reionized if Chii is lower than the red line, 
if a compact morphology is used (template 5), or lower 
than the blue line if an extended one is considered (tem- 
plate 4). The puv is computed integrating the LF down 
to M1500 = — 10 magnitude. The solid line at Chii = 1 
indicates a homogeneous IGM, and we must consider only 
the Chii > 1 region, since galaxies at z ~ 7 are formed 
in biased density regions of the Universe. The dashed line 
shows the value Chii = 3. In practice, only a combination 
of (fescCnn) parameters inside the shaded regions corre- 
sponds to a z ~ 7 IGM kept ionized by LBGs only. 



(G2.1408 of lCastellano et al. 2010a|) . If the faint end of the 
LF at z=7 is not so steep, as we have inferred here for the 
compact morphologies, the contribution of fainter galaxies 
(-19 < Mi 5 oo < -10) is only 2.4 times that of Afi 50 o < 
— 19 LBGs. Thus, the present uncertainty on the steepness 
of the faint end of the LF translates in a large uncertainty, 
of a factor 2-15, on the best fit of the luminosity density, 
SFR density, and reionization rate at z=7. 

Following Wind horst et al 20021 it is plausible that 
fainter galaxies have smaller physical sizes: in this case 
a combination of high escape fraction (f esc > 0.2) and 
small dumpiness {Chii < 3) is required for the reion- 
ization of the Universe. It is worth noticing that the 
recent estimates of f esc for L > L* LBGs at < 
z < 3 are < 15 (|Bridge et al. 20101 ICowie et al. 20101 
ISiana et al 20101 IVanzella et al. 2010[) . implying a consis- 
tent evolution of the galaxy escape fraction going to fainter 
luminosities or to higher redshifts. 

Moreover, the knowledge of the reionization processes in 
the high-z Universe is currently far from clear. For example, 
early theoretical works (i.e. IGnedin fc O strikcr 1997) have 
found high values for the dumpiness, Chii ~ 10 — 30. More 
recently, IPawlik et al. 20101 studied the effect of photoion- 
ization heating by a uniform ultraviolet background and 



found that photo-heating can strongly reduce the clumping 
factor of the IGM because the increased pressure support 
smoothes out small-scale density fluctuations. The common 
expectation is that photoionization heating should provide 
a negative feedback on the re-ionization of the IGM be- 
cause it suppresses the cosmic SFR by boiling the gas out 
of low-mass DM halos. However, since it also contribute to 
the reduction of the clumping factor, the net result is that 
it is easier to keep the IGM ionized. Photo-heating there- 
fore also provides a positive feedback for the reionization 
processes. Neglecting the effects of photoionization heating, 
Pawlik ct al. 2010 predict a dumpiness of the IGM of the 
order of 7-8 at z=7, and in this case our limit to the LF rule 
out a completely ionized IGM even for f esc — 1 if all the 
galaxies at z — 7 resemble our morphological template 5 
(compact galaxies), while an f esc ~ 0.3 is still sufficient 
assuming template 4 (extended galaxies). If instead the 
photoionization heating is considered, a small dumpiness 
Chii < 3 is expected, in agreement with measurements by 
Bolton & Haehnelt (2007). In this case, an escape fraction 
of ^50% is enough to keep the Universe re-ionized, even in 
the most pessimistic case (of compact morphology) . 

6.2. The reliability of our LF estimate 

The uncertainties on the estimate of the z = 7 LF are re- 
lated to the method used to derive it. Errors can be due 
to the limited depth of the images used, the color criteria 
adopted, the presence of possible interlopers (cool dwarf 
stars, low-z red galaxies, transient objects), spurious detec- 
tions, especially when going faint. All these features have 
been discussed in the previous sections, showing that de- 
tailed simulations are able to take into account the effects 
of different color criteria, and that a precise color cut can be 
identified with the aim of separating the few high-z galaxies 
from the much frequent interlopers. Moreover, the details 
on the steepness of the faint end of the LF are strongly 
influenced by the morphological templates adopted for the 
simulations. 

There are however two issues that cannot be corrected 
with the simulations only: the presence of transient objects 
within the z-dropout candidates and the role of the cosmic 
variance when dealing with deep, pencil beam surveys. In 
the near future e.g., the HUDF field will be re-observed 
with WFC3 in order to go deeper than the present obser- 
vations, with 34 additional orbits splitted in three filters 
(Y, J, and H), complementing the existing 60 orbits dis- 
cussed here. With these data, it will be possible to assem- 
ble a z-dropout sample clean from variable contaminants. 
The discovery of a clear transient source in the HUDF data 
(see FigJT|) and two possible transients in the Hawk-I data 
by ICastellano et al. 2010al is more than an order of magni- 
tude higher than the 0.06 SNe per WFC3 field estimated 
in lOesch et al. 20101 It is possible that some of the selected 
z=7 candidates will be rejected as interlopers when a sec- 
ond epoch survey will be executed. 

Surveys are usually bounded by a limited amount of 
telescope time, and thus they have to compute a trade 
off between the covered area and the imaging depth. This 
has deep implications on the number of objects recovered, 
on their absolute magnitudes and the sampled Universe 
volume. Simple observational estimates of galaxy number 
counts in finite volumes are thus subject to the uncertain- 
ties due to cosmic variance, arising from underlying large- 
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scale density fluctuations. This may be particularly true at 
high-z, where galaxies are expected to be much more clus- 
tered than at lower-z. 

We have used the Millennium simulation, using the real- 
ization of Kitzbichlcr White 20071 to investigate the im- 
pact of the cosmic variance on the number of galaxies re- 
covered in the two WFC3 fields analysed here. We extract 
randomly from the simulation of Kitzbichlcr & White 20071 
10 4 pencil beam realizations, recording for each pointing 
the number count N of galaxies at 6.5 < z < 7.5 down 
to Misoo = -17.5 and -19.3 for the HUDF and ERS, 
respectively. The cosmic variance is computed as &q V — 
°m — 1 / < N >, where the first term is the measured vari- 
ance and the second term is the poissonian contribution. We 
found that the expected uncertainty due to cosmic variance 
is <7 C v/ < N >=~ 34% and - 38% for the HUDF and ERS 
field, respectively. Combining the two fields, the resulting 
cosmic variance yields ~ 28%. The se results are in agree- 
ment with the predictions based on ISomerville et al~2 004, 
extrapolated to z ~ 7. 

As shown in IRobcrtso n 20101 since the errors on M* , 
$*, and a are correlated, this does not translate directly 
in an additional ~ 30% uncertainties in the LF parame- 
ters. However, the uncertainties on the luminosity density 
or on the UV ionizing photons depends mainly on the total 
number of galaxies, so the ±30% additional uncertainties 
can be applied to these quantities to take into account the 
effects of cosmic variance. 



7. Summary 

We have analysed two WFC3 fields, the HUDF and ERS 
surveys. More than 40 bona fide z-dropout galaxy candi- 
dates have been selected on these two areas, down to an 
observed magnitude of J - 29 (AB). The depth of ACS 
images observed with bluer filters (BVI) ensures a robust 
and well controlled selection of the high-z galaxy sample, 
cleaning the candidates from the large number of interlop- 
ers (i.e. low-z early galaxies, dusty starbursts or cool stars). 

We have explored different color criteria for the selec- 
tion of z ~ 7 galaxies, picking up three methods from the 
wealth of methods proposed recently in the literature. The 
numbers of the resulting candidates are very different, with 
13, 6, and 20 candidates following the A, B, and C criteria 
(described in Sect. 3) in the HUDF field, while in the ERS 
we find 14, 6, and 22 z-dropout candidates, respectively. 
We have shown here that, applying these different methods 
both to the observed surveys and to the simulated samples, 
the resulting LFs are similar. 

On the other hand, the output LF depends critically 
on the morphological templates used to simulate the com- 
pleteness of the observed candidates. In Section 4 we have 
proposed 10 different options for a reliable simulation, span- 
ning from point like objects to very extended galaxies, both 
observed and synthetic. The resulting LF is extremely sensi- 
tive to the adopted morphological templates, especially the 
parameter a governing the steepness of the faint end. We 
have shown that the steep LF found bv lBouwens et al~ 2011 
is possibly due to the choice of a set of templates with typ- 
ical half light radius larger than the one of the observed 
candidates at z=7. If a more co mp act template is adop ted 
(as used in IMcLure et al. 20091 and IWilkins et al. 2011[) , a 
milder slope is derived. 



The resulting best fit to the z ~ 7 LF is Log(<&*) = 
— 3.13±g;H, M* = — 20.14±g;f2, and a = -1M±° 7 5 1. This 
implies a luminosity density (down to an absolute magni- 
tude of -10) of p uv (-10) = 9.2tfi 8 • 10 25 erg/s/Hz/Mpc 3 
and a star formation rate density of SFRD(— 10) = 
0.017i ; 010 M g /yr/Mpc 3 , if a compact morphology is used 
(template 5). If a more extended morphology is adopted, 
the resulting best fit to the z ~ 7 LF is Log(&*) = 
— 3.22±g;||, M* = -20.3Q±g;|^ and a = -2.02±gj£. This 
implies a luminosity density (down to an absolute magni- 
tude of -10) of puv(-lO) = 4.5tfe ■ l0 26 erg/s/Hz/Mpc 3 
and a star formation rate density of SFRD(— 10) = 
O.OSStlf^Mo/yr/Mpc 3 . 

The contribution of z=7 LBGs to the hydrogen ionizing 
photon density Ni on depends on the last two parameters, 
and the threshold to keep the Universe ionized at z ~ 7 
depends also on f esc , the escape fraction of ionizing pho- 
tons, and on the dumpiness of the IGM, Chii- Limiting 
the dumpiness to Chii > 1, we thus have f esc > 0.17 in 
order to have the Universe reionized by z=7. Inferring a 
small clumping factor for the IGM at high-z (Chii < 3 
at z > 6, Bolton & Haehnelt 2007), an escape fraction of 
fesc ~ 0.5 is enough for the reionization of the IGM (but 
not for the more clumpy galaxy formation regions) , assum- 
ing a LF derived from template 5 (compact morphology). 
If the IGM was reheated at z > 9, the observed population 
of LBGs at z=7 might be enough to keep the IGM ionized, 
given that the escape fraction is larger than 50%, since a 
dumpiness of Chii ~ 2 — 3 is expected, according to recent 
theoretical predictions (jPawlik et al. 2010]) . 

In the near future, two large surveys will provide 
further information and better statistics to the z ~ 
7 related studies: the ground based wide NIR survey 
ULTRAVISTA0 (Ultra Deep Survey with the VISTA tele- 
scope) on the COSMOS field and the deep space based 
survey CANDEL£|f| (Cosmic Assembly Near Infra-red Deep 
Extragalactic Legacy Survey). In particular, the CANDELS 
survey will have a deep impact on the z=7 study since it will 
be able to clean out the present high-z sample from vari- 
able objects, to reduce significantly the cosmic variance, to 
sample in detail the break of the LF at Afisoo ~ —20, and 
to provide a detailed distribution of the half light radii of 
the candidates, strongly needed to derive a reliable incom- 
pleteness function from the simulations. 

In the meanwhile, spectroscopic confirmations of 
z=7 (and beyond) candidates wi th optical/NIR in- 
strumentations from the ground (Fontana et al. 2010, 
IVanzella et al. 201 1|) . with JWST from space or through 
sub-mm emission lines with ALMA will open new frontiers 
in the study of the first galaxies. 
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